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Abstract: We approach the following question: if supersymmetry is discovered, 
how can we select among different supersymmetric extensions of the Standard Model? 
In particular, we perform an analysis of the sparticle spectrum in low-energy string 
effective theories, asking which observables best distinguish various scenarios. We 
examine scenarios differing by the fundamental string scale and concentrate on GUT 
and intermediate scale models. We scan over four parameters (two goldstino angles, 
tan j3 and the gravitino mass) in each scenario, finding ratios of sparticle masses 
that provide the maximum discrimination between them. The necessary accuracy 
for discrimination is determined in each case. We find that the required accuracy 
on various sparticle mass ratios is at the few percent level, a precision that may 
be achieved in future linear colliders. We also map out phenomenologically viable 
regions of parameter space. 
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1. Introduction 

There is great expectation in the high energy physics community of the possibility 
of discovering evidence of low-energy supersymmetry (SUSY) during the next few 
years. The smoking-gun signatures of production and detection of super-partners 
could be observable at the Tevatron, the Large Hadron Collider and a future linear 
collider facility. If such signatures are detected, we may enter a new era of high energy 
physics with closer contact between fundamental theory and experiment. Since there 
are many supersymmetric models, the issue may then turn from discovering SUSY 
to selecting and eliminating the different supersymmetric models that have been 
proposed over the years. It is therefore useful to examine ways of comparing the 
different models. Experiments may start to give us information not only about 
particle content but also about the mechanism of SUSY breaking and (eventually) 
the messenger of this breaking. 

Many argue that the best motivated models of low energy SUSY are those that 
can be derived from string theory. Even though there is a plethora of possible specific 
string models, we can attempt to perform a fairly model independent analysis using 
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some general string scenarios. There are several ways to parameterize these scenarios. 
In recent years it was realised that the underlying scale of a fundamental string 
theory can be different from the Planck scale of 10 19 GeV if the observable fields are 
confined to a brane within a higher dimensional world. Since the size of the extra 
dimensions is not necessarily fixed in such a theory, we have the freedom to argue 
for different values of the string scale. Two such possibilities with several indications 
in their favour are the GUT scale Mqut ~ O(10 16 ) GeV |Ij] and the intermediate 
scale of around Mi = 10 11 GeV [fj, [3|. The GUT scale is favoured by the data in 
the MSSM-desert gauge unification picture 0. The intermediate scale is motivated 
by a natural solution to the strong CP problem, the scale of neutrino masses || 
and a natural supersymmetry breaking scale in gravity mediated supersymmetry 
breaking scenarios. But in this case, assuming an MSSM-like low energy spectrum, 
the gauge couplings evolved from the data at the electroweak scale M z to Mj do 
not meet, contrary to the naive string prediction. One must address the issue of 
gauge unification in the intermediate scale models, and two possibilities have been 
identified: one possibility is to achieve precocious unification || through the inclusion 
of additional matter fields on top of the MSSM content. Explicit intermediate scale 
string models have sometimes exhibited the existence of such superfields extra to 
the MSSM [|J. Another possibility is mirage gauge coupling unification 0. In 
this picture, gauge unification at the GUT scale is merely an illusion created by 
string loop effects; the field theoretic gauge couplings are set non-universally at the 
(intermediate) string scale. Here, no additional matter fields beyond the MSSM are 
necessary. 

We can imagine having experimental data on supersymmetric particles and would 
like to see if they may be able to determine the favoured string scale, or other im- 
portant parameters of the fundamental theory. In this way we would start to obtain 
relevant information about the fundamental theory at a large scale from measure- 
ments at low energies. Motivated by this possibility we ask to what extent we can 
discriminate scenarios for the underlying theory. See || for related work. 

Supersymmetric phenomenology is extremely complicated. The sparticles un- 
dergo cascade decays which provide various signatures in experiments. Such signa- 
tures notoriously depend not only on the model of SUSY breaking, but also upon its 
parameters. If a SUSY breaking parameter changes, mass differences in the sparti- 
cle spectrum can change sign. Decay channels in the cascade sparticle decays then 
switch on and off, and the identified signature is easily lost. Such is the complexity 
of supersymmetric phenomenology that experimental studies of future colliders have 
typically focused on a few points in parameter space |7j . 

To avoid getting bogged down in such complexity and model dependence, and 
in order to start the phenomenology ball rolling, we will assume that some of the 
superparticles' masses have been measured. We will try to gain information on the 
high energy theory from their values. Explicitly, we will examine three different 
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scenarios: 



• String scale at the GUT scale Mqut ~ O(10 16 ) GeV, defined by the scale of 
electroweak gauge unification gi^Mcur) = 92(Mgut)- 

• Intermediate string scale (Mj = 10 11 GeV) with extra leptons to achieve gauge 
coupling unification at Mj, which we will refer to as early unification. 

• Intermediate string scale (Mj = 10 11 GeV) with mirage unification |J. The 
particle spectrum is assumed to be as in the MSSM, and the gauge couplings 
at Mi are unequal. 

To predict sparticle masses from these scenarios, we must solve the renormalization 
group equations starting from a theoretical boundary condition parameterized by 
the string scale, the goldstino angles, tan/? and the gravitino mass 7/13/2- Constraints 
from experiments and cosmology (if a version of i?-parity is conserved, as assumed 
here) restrict the models. 

Previous studies have been performed on intermediate string-type scenarios. Two 
of the authors || previously mapped out the spectrum and naturalness parameter for 
the dilaton dominated scenario in the GUT, early unification and mirage scenarios. 
Baek et al extended this study to include bounds from B — ► X^, dark matter and g — 
2 of the muon || . The early unification spectrum and relic neutralino density has also 
been investigated [H| when one of the goldstino angles varies. However, the effects 



of moduli fields on the SUSY breaking was neglected in this study. Dark matter 
observables were also considered for mirage or early unification scenarios in ref. 
Dilaton domination (corresponding to a fixed limit of one of the goldstino angles) is 
in violation of a charge and colour breaking bound for GUT scale unification [[1^1 , 
but it was found that this bound does not restrict the intermediate scale case |§. 

Here we extend previous investigations in three directions. Firstly, we scan over 
both of the goldstino angles. Secondly, we look for the combinations of sparticle mass 
measurements which provide the most discrimination between the string scenarios, 
once parameters are scanned over. Thirdly, we approach the question: is it possible 
to select a string scenario from a knowledge of the masses of the different sparticles? 
We re-phrase this question as: what accuracy in sparticle mass measurements is 
required to select a string scenario? 



2. Soft SUSY Breaking Terms 

The Type I string models considered here are based on orientifold compactifications 
of type IIB strings, and share a number of similarities with the well studied phe- 
nomenological models derived from compactified heterotic string theories |13|]. Both 
contain a dilaton superfield S, and moduli fields Tj connected with the size and shape 
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of the extra dimensions. These two fields contribute to SUSY breaking when their 
auxiliary fields acquire vacuum expectation values (VEVs). We consider one overall 
modulus T and later one blowing-up mode M as a parameterization of the stringy 
SUSY breaking. The respective contribution of F s and F T to the SUSY soft break- 
ing terms can be parameterized in a goldstino angle 9 (where sin 9 = 1 corresponds 
to dilaton domination and cos 9 = 1 denotes moduli domination). 

One difference between heterotic and type I string models is that the string scale 
is not constrained to be near the Planck mass. Furthermore, the set of moduli fields 
connected to the blowing-up modes, M a , play a more relevant role for SUSY breaking 
in Type I models than in the heterotic ones. They contribute explicitly to the gauge 
kinetic coupling and therefore their F-terms may induce gaugino masses and the 
other soft breaking terms. M mixes in the Kahler potential with the modulus field 
T. Accordingly it is convenient to introduce a second, (F T , F M ) mixing angle in 
analogy to the usual dilaton-moduli mixing angle 9. To summarise, 
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where F total = V f s + f t + f m- 

The introduction of 9 and provides a convenient way to parameterize the 
influence which the various VEVs have on the soft terms. We will generalise the 
results of ref. || where only the dilaton domination case was considered extensively. 
Here, we examine a wider 9 and range but leave out 9 < 10°. For such small values 
of 9, the string-induced high scale boundary conditions on the soft breaking terms 
reach a similar magnitude as the one-loop order anomaly mediated SUSY breaking 
terms, which therefore cannot be neglected. The complete set of soft breaking terms 
of combined anomaly and gravity mediation has not yet been computed fl4" |. We 



therefore leave this slice of parameter space to the future, when such terms might 
be included in an analysis. Like the analysis in ||, we assume that the full standard 
model gauge group arises from a single brane and that SUSY breaking is dominated 
by the F terms of the S, T and M fields. This results in the following soft breaking 
terms in the T + T —>■ oo limit: 

2.1 Scalar masses 

All scalar masses receive the universal soft SUSY breaking term 

ml = Vo + m 2 3/2 (l - \C 2 cos 2 9 sin 2 0^ + O ( (T ^ )2 ) , (2-2) 



where C = A /l + 3 J^? , and Vq is the vacuum energy, k depends on the form of the 
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Kahler potential and can be set constant 1 ||. We set k = 3 in this article, which 
avoids negative scalar mass squared values at the string scale. Higher values of k 
are possible, this would lead to a weaker dependence of the scalar soft masses on 
the goldstino angles. As usual we will take C = 1 corresponding to V = 0, i.e. a 
vanishing cosmo logical constant. 

2.2 Gaugino masses 

The gaugino masses are equal to 

M a = V3Cm 3/2 ^- (sin 9 - ^a GUT cos 9 sin + cos ^ ) +0 (<— ^j) ■ 

(2.3) 

Here 5gs is the Green-Schwarz term coming from anomaly cancellation (like in com- 
pactified heterotic string models). Its value is a model dependent negative integer 
of order 0(— 10). We will fix its value to —10 from here on 2 . Note that, in general, 
the gaugino masses have a non-universal boundary condition at the string scale. If 
the string scale is intermediate, a a =x,2,z 7^ &gut provides non-universality in equa- 
tion (|2.3| ) from the first term. The second term provides non-universality at the 
string-scale through one-loop stringy effects away from the dilaton dominated limit 
(9 = 90°). This term is proportional to the model-dependent parameters s a . In order 
to make mirage unification possible, s a = 2n/3 a was chosen where (3 a are the usual 
MSSM renormalization /5-function coefficients (33/5, 1, —3). 

2.3 A-terms 

Under the assumption that the Yukawa couplings are moduli- independent, the tri- 
linear couplings are 

A aPl = -V3Cm 3/2 ( sin 9 + cos 9 cos <f>K'^ + O (^^j . (2.4) 

Here K = K(M + M* - 5 GS \og(T + T*)), the M- and T-dependent part of the 
Kahler potential. One can set K =0 M, since all fields are assumed to be in the 
minimum of the potential, for which the argument of K vanishes. 



1 Usually in heterotic models the scalar masses are non- universal, since the coefficient of cos 2 9 
depends on the so-called modular weights of the corresponding matter fields |f5|, |i~3f . If we con- 
sidered several T fields here, we could also expect a non-universal behaviour, parameterized by 
different values of the constant k. However, for a single T field we get the universal behaviour 
shown in ( |2.2| ) , which amounts to saying that in these models all matter fields have modular weight 
— f . For a recent discussion of non- universal soft terms in a different context, see for instance Jf6| . 

2 Different values for <5qs affect the dependence of M a on the goldstino angle 4>. While this is 
negligible for \Sgs\ < - 1 - where the minimum of M a stays very close to 9 = on the (9, <f>) plane, 
higher values of <5gs lead to a significant perturbation away from 9 = 0. For a discussion of variable 
S GS , see pi. 
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3. Renormalization Group Analysis 



In order to analyse the sparticle spectra we use the ISASUGRA part of the ISA JET 7.51 
package |18| to obtain the SUSY spectrum starting from the high-energy boundary 



conditions detailed above. 
3.1 Constraints 

We use the following experimental constraints to limit the scenarios 



m^o > 45 GeV m~± > 103 GeV m ho > 113.5 GeV. (3.1) 

We will also use constraints from the recently measured muon anomalous mag- 
netic moment = (g — 2)/2. The experiment E821 at Brookhaven National Labo- 
ratory (BNL) reported the measurement 



a M = (11659202 ± 14 ±6) x l(T iU . (3.2) 

The precision in the measurement already significantly constrains the SUSY param- 
eter space and will improve in the future. We will constrain all investigated models 
to be within 90% CL of the central value of the measurement 22]] , i.e. 



-4.2 x 10~ iU < 5a„ < 41.3 x 10~ iU . (3.3) 

The final constraint is the absence of a charged lightest supersymmetric particle 
(CLSP). This constraint must be applied when the LSP is stable on cosmological 
time scales. This is the case here because we implicitly assume that a version of 
.R-parity holds, as recently found in specific string models |23[ . 



3.2 String scale boundary conditions 

As mentioned above, we consider three overall scenarios: GUT scale unification (with 
the string scale set to Ms = 2 x 10 16 GeV), early unification and mirage unification 
(both with Ms = Mi = 1 x 10 11 GeV). The MSSM spectrum is assumed in each 
case, except for early unification where 2 x Ll + 3 x Er vector-like representations are 
added to the MSSM in order to achieve gauge unification at Mi M. It is assumed 
that these extra representations have negligible Yukawa couplings. We add their 
effect to the MSSM gauge /5-functions above 1 TeV. 

The string-scale boundary conditions on the soft terms were obtained from equa- 
tions ( p.2[ - p.4|) , with m 3 / 2 , 9 and 0, as well as tan/3 as free parameters. Throughout 



the whole analysis /z > and m t = 175 GeV were assumed |19| . Negative /i leads to 
a negative Sa^, which is limited from equation [O] to be small in magnitude. This 
means that, for a given value of tan/?, the sparticles must be heavy in order to sup- 
press their contribution to 5a^ (compare figs. |3] and |j). In this limit, effects of the 
sign of [i upon the mass spectrum are suppressed. We can therefore safely ignore the 



- 6 - 





9 


<P 


™3/2 


tan (3 




m t 


V 


c 


fc 


<5gs 


M GUT 


Mi 


range 


10-90° 


0-90° 


50-1500 


2-50 


> 


175 





1 


3 


-10 


2xl0 16 


10 11 



Table 1: Summary of parameters. The first four parameters are scanned over, and their 
range is detailed. For the others, their value is kept constant except for which is con- 
strained to give the correct value of Mz- All massive parameters (m 3 /2,mt, Mi, Mqut) 
are given in units of GeV. 



Ii < case because its resulting spectra will be included in our /j > results. This 
also means, though, that the sign of /J cannot be determined from these spectra, 
unless the sparticles are too light to be compatible with /j < for any tan f3. 

In the GUT case, a a = ogut = 1/25 was used in equation (|2.3|) for all gauge 
groups a = 1, 2, 3. In the mirage unification case, however, the gauge couplings only 
appear to unify at the GUT scale. In reality they are set at the intermediate scale 
Mi = 10 11 GeV to non-universal values 3 . In our analysis, 

111 

aAMA = , a 2 (Mi) = , a 3 (Af 7 ) = (3.4) 

' 37.6' V ; 27.0' V ; 19.8 V ; 

is the set of values which remains reasonably stable under iterations of inserting them 
into equation ( |2.3j ) and running the RGEs again. 

For each of the three scenarios that we consider, we scan over the free parameters 
9, </>, tan/? and m 3 / 2 . 9 = 90° corresponds to the dilaton domination case considered 
in M, and the SUSY soft mass terms are independent of </> in this case. As men- 
tioned above, values of 9 < 10° were not considered, since in this case tree level soft 
masses would vanish and it would not be possible to neglect the one-loop effects from 
anomaly mediated SUSY breaking, which might modify the soft terms considerably. 
Table [I] shows the default ranges and increments of scanned parameters. It also sum- 
marises the other parameters that were kept constant. When different values were 
used for a plot, we will mention it explicitly. 

3.3 Sparticle spectra 

In order to get some feeling for the effect of the goldstino angles, we first illustrate 
their effect upon the sparticle spectra. We choose the mirage unification scenario as 
a first example. The other scenarios do differ from these spectra quantitatively, but 
the effect of the goldstino angles is similar in each case. In figure [1], we show the 
variation of the mirage unification spectrum with 9 for tan/5 = 30, 7713/2 = 300 GeV 
and = 60°. As mentioned above, 9 = 90° corresponds to the dilaton dominated 
limit. From the figure, we see that larger 9 increases the splittings between sparticles. 
Larger 9 corresponds to larger gaugino masses in eq. (|2.3| ) and the larger gluino mass 

3 We note that ref. Q assumes gaugino universality, which leads to different conclusions for the 
mirage scenario, especially for the charged LSP constraints. 
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Figure 1: Variation of mirage unification sparticle mass spectrum with 9. We use eft = 60°, 
tan/3 = 30 and m 3 / 2 = 300 GeV. The key on the right-hand side of the figure details the 
flavour of sparticle. 



raises the other coloured sparticle masses in the running from Mi to Mz- The 
weak gauginos also show a milder change from the larger values of M lj2 . Sleptons 
are largely unaffected by the variation of 9, demonstrating the fact that the mild 
increase in mo in eq. (|2.2|) from increasing 9 is a relatively minor effect. 

The ordering of superparticle masses is largely unaffected by changes in 9, except 
for the lightest chargino, which crosses the e^, f 2 and ej, lines 4 and the lightest stop, 
which crosses the heavier chargino line. Each ordering can correspond to different 
cascade decay channels, and might be used to restrict 9 once these decays are ob- 
served. For example, when the chargino crosses the cr line, nothing much changes 
because the lightest chargino does not have significant couplings to right-handed se- 
lectrons. However, decays into neutrinos and or f 2 (which can have a significant 
left-handed component for large tan/3) become viable for 9 > 70°. 

In figure |2]a, we show the variation of the mirage unification sparticle spectrum 
with for tan/3 = 30, m 3 / 2 = 300 GeV and 9 = 30°. Here, increasing has a 
smaller effect than 9, but generally decreases the splittings between sparticle masses. 
The effect is relatively small because the gaugino masses are only sensitive to 
through a suppression factor in eq. dO[) . Most of the sparticle masses decrease 
with increasing 0, except for the lightest neutralino (which becomes heavier) and 
the lightest chargino (which remains roughly constant). This can be understood by 



*It is difficult to discern in the figure that the fi and are quasi-degenerate. 
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Figure 2: Variation of sparticle mass spectrum with <p for (a) mirage and (b) GUT-scale 



unification scenarios. We use 9 = 30°, tan/3 = 30 and m 



3/2 



300 GeV. The key on the 



right-hand side of the figure details the flavour of sparticle. 



considering the effect of increasing <p on the scalar masses in eq. (|2.2| ). The figure 
shows that the model would be ruled out for > 70° because of a stau LSP. Aside 
from this, important changes in the ordering of sparticle masses occur between the 
lightest stau and chargino and the heavier chargino and the lightest stop. The heavier 
stop becomes heavier than the gluino for > 50°. This is potentially important 
because future e + e~ machines can produce and measure gluinos more easily through 
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the decays of pair-produced squarks. In this example, however, the mass difference 
between gluino and the heavier stop is not sufficient to produce an on-shell top quark 
in a decay. Thus such a process would be highly phase-space suppressed. 

Figure shows the variation of the GUT scale unification sparticle spectrum 
with for tan/? = 30, 7713/2 = 300 GeV and 9 = 30°. While quantitative differences 
exist between figures F2]a and FJb, it is clear that the two spectra are qualitatively 
very similar. There is, however, some re-ordering of sparticle masses between the 
GUT-scale and mirage-scale scenarios for certain values of 0. For example, at = 
20°, the heavy chargino and light stop are interchanged. It would be difficult to 
discriminate the two scenarios on the basis of decays implying a certain ordering of 
the lightest stop and chargino since there is no a priori information about 0, although 
it might be possible to fit the masses to certain values of 9, if the experimental 
accuracy were good enough. For larger values of 9, the variation of the spectrum 
with significantly decreases because the terms that depend upon it in the boundary 
conditions eqs. (^.2|)-(PT^D become smaller and sub-dominant to the other 9 dependent 
terms. 

3.4 Maps of parameter space 

We now explore the parameter space, including constraints from the sparticle mass 
limits in eq. (|3.1| ), requiring that the g—2 anomalous magnetic moment be within 90% 
CL of the central measured value, requiring correct electroweak symmetry breaking 
and imposing a neutral LSP. In figure |3|, we present viable regions of tan f3-m 3 / 2 
parameter space for various points in the goldstino angle space and /i > 0. Dilaton 
domination has been well covered in the literature || || and so we present other 
more general regions of the parameter space alongside the dilaton dominated limit 
in figure ||b. We see from the variation in the viable regions from figures [3|a-|3|d that 
they are dependent upon the goldstino angles, which vary between each of the figures. 
For 9 = 30° and = 30°, we see from figure |3]a that there is little dependence of 
the viable region upon the scenario assumed. The other three figures show a larger 
difference between the three scenarios, apart from for = 60°, 9 = 30°, which does 
not distinguish the GUT and mirage unification regions significantly, as can be seen 
from fig. |3|c. In figure [|, we show the available parameter space for fi < 0. By 
comparing figs. [3] and |4], we see that < admits less parameter space. This is 
because the g — 2 constraint restricts sparticles to be heavy for /i < in order to give 
a small contribution to the magnitude of the anomalous moment. The lower bounds 
upon sparticle masses are much stronger for high tan (3, where the contributions to 
the anomalous moment are largest. Figure (| illustrates that (for low tan/3), there is 
a lower bound m 3 / 2 > 400 GeV. 

There are often also bounds upon tan (3. Low tan /3 is ruled out by the direct lower 
bound upon the Higgs mass from LEP2 |fT9"| , whereas high tan/5 is either ruled out 
by the constraint of a neutral LSP, or of not too large a supersymmetric contribution 
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Figure 3: fi > maps of parameter space for the early unification (red crosses), mirage 
(blue dots) and GUT unification scale (green diamonds) scenarios. Here, we plot viable 
regions of tan parameter space for fixed values of the goldstino angles, after all 

exclusion limits were applied: (a) 6 = 30°, <j) = 30°, (b) 6 = 90°, <fi =arbitrary, (c) = 30°, 
(j) = 60°, (d) 9 = 60°, <f> = 30°. 



to g — 2 of the muon. For instance in the dilaton domination fi > (fi < 0) case we 
can see that the GUT scale scenario restricts tan/3 < 40 (44), the early unification 
scenario requires tan/3 < 30 (28) and mirage unification constrains tan/3 < 46. For 
the lower value of 9 = 30° in the [i > case, the tan j3 bound disappears for = 30° 
in all scenarios as in figure ^|a, but reappears for higher </> = 60° in the early and 
GUT-scale unification scenarios (figure |3]d). 

In figures [||a and [|b, we show the viable regions in goldstino angle space for 
fixed m 3 / 2 = 400 GeV and tan/3 = 20,40 respectively. For lower tan/3 = 20, there 
are viable dilaton-dominated cases for each scenario. However, for tan/3 = 40 in 
figure |5|b, dilaton domination (6 = 90°) is not viable except for the mirage unification 
case. The bounds <p < 50° an d < 70°, 60° apply to the GUT-scale and early 
unification scenarios respectively. For = 0°, there is still a large contribution to 
unification-scale sparticle masses from terms involving 6, so a lower bound on <p does 
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Figure 4: < maps of parameter space for the early unification (red crosses), mirage 
(blue dots) and GUT unification scale (green diamonds) scenarios. Here, we plot viable 
regions of tan (3-m^/2 parameter space for for fixed values of the goldstino angles, after all 
exclusion limits were applied: (a) 9 = 30°, <\> = 30°, (b) 9 = 90°, (ft =arbitrary, (c) 9 = 30°, 
(j) = 60°, (d) 9 = 60°, <p = 30°. 



not exist for general 6 from sparticle mass bounds. 
3.5 Discriminating ratios 

We have searched through ratios of masses of sparticles to see which ones provide 
the largest discrimination between the three scenarios examined. The ideal situation 
would be to find ratios that are predicted to be completely disjoint regions for the 
three scenarios. Then, if the mass ratios were determined in an experiment, a com- 
parison with these disjoint regions should reject or confirm one scenario. In general, 
this may not be possible so we resort to finding ratios which discriminate between 
some of the scenarios in some cases. Where disjoint regions are found, we will focus 
on the accuracy required in any mass ratio measurement that would be enough to 
discriminate two scenarios. This accuracy we will define as follows: supposing the 
central value of a measurement of a ratio lies on the edge of one of the regions. The 
fractional accuracy in each direction that would exclude each of the other scenarios 
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Figure 5: Maps of parameter space for the early unification (red crosses), mirage (blue 
dots) and GUT unification scale (green diamonds) scenarios. Here, we plot viable regions 
of 9 -cj) parameter space for fixed values of and tan (3: (a) 7713/2 = 400 GeV, tan (3 = 20, 
(b) m 3/2 = 400 GeV, tan (5 = 40 



is taken. Then, we examine all points lying within a region and quote the lowest 
fractional accuracy on each ratio required to separate the regions. One ratio on its 
own does not provide enough separation once the parameters in each scenario are 
scanned over; we therefore resort to considering combinations of two ratios in order 
to see if the regions of each scenario are separated in two dimensional ratio space. 

There are a few combinations of ratios that provide near separation. We display 
the best two such combinations in figure || For each point in the parameter scan 
described in table |l], that satisfies the constraints in eq. (3.1) as well as the neutral 



LSP constraint, one point appears on the plot. Here, we have taken input parameters 
at random within the ranges defined in table |l] for 10000 points (as is the case for 
the rest of this subsection). Each of the three scenarios were scanned over and 
shown on the plots. We have separated figure |6| into two cases: either the squark 
flavours/handedness are known, as assumed in figure ^|a, or a squark mass averaged 
over the first two generations and handedness is used, 



1 



in 



*»» - g ( m u L + m u R + m c L + ™~c R + m St + m SR + m~ SL + m- SR ) , (3.5) 

as in figure |6|b. For our numerical analysis we assumed degeneracy between the first 
two families, and therefore in all our figures, and m^ L can also stand for mg £ 
and rris L respectively. The reason for showing both figures is the differing ability of 
hadron and lepton collider machines. It will be impossible for a hadron machine like 
the LHC to distinguish between the flavour of the almost-degenerate first two families 
of squarks. It is also difficult to see how a mass measurement of different handedness 
of squarks might occur. The third family typically are split from these in mass, and 
have additional heavy flavour tagging identification for their decay products and so 
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(b) 




Figure 6: Highly discriminating ratios of sparticle masses. Each plotted point corresponds 
to one data point in four-dimensional parameter space (9, (j), tan/3, m.3/2) that is not 
excluded by the imposed limits. In (b), rriq avg is as defined in the text. Pluses correspond 
to the early unification scenario, crosses to GUT-scale unification and stars to mirage 
unification. 



it should be possible to separate these. On the other hand, possible future linear 
e + e~ colliders are expected [23| to provide very precise measurements of mg L and 
mg R . Because a future linear collider could scan to an energy that is sufficient to 
produce a particular mass of squark, it will be possible to separate up-type from 
down-type squarks provided there is enough energy to produce them, and identify 
which type squarks are produced from their cross-section. Handedness could either 
be determined by examining the decay chains, or by using polarised beams to favour 
production of one handedness or another. It is conceivable that squarks will only be 
produced at the LHC, and are too heavy for a particular linear collider facility to 
produce. It is for this reason that we provide discrimination plots for both individual 
and averaged squark masses. The importance of possessing individual squark mass 
measurements can then be assessed. 

In figure ^|a there is a large overlap between the early unification and mirage 
scenarios' predictions in the range 0.75 < m<j/m Si < 1.15, but outside of this range 
separation is possible. GUT-scale unification can usually be discriminated from the 
other two scenario ranges, provided we measure the mass ratios with 1% accuracy. 
In Fig. Pd, we see that the effect of averaging over the squark mass measurements is 
to significantly decrease the discrimination between the mirage and GUT unification 
scenarios. There is significant overlap in the allowed ranges for all three scenarios, 
but that there are also regions of predicted mass ratios which can are only consistent 
with one particular scenario towards the bottom-left or top-right of the plot. If the 
empirical uncertainties on mass ratios are small enough and the underlying model 
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Figure 7: Highly discriminating ratios of sparticle masses for dilaton domination (6 = 
90°). Each plotted point corresponds to one (m 3 / 2 , tan/3)-pair in the scan not excluded by 
any limits. Green crosses are predicted by the GUT unification scale, early unification is 
denoted by red crosses, and blue stars are valid for the mirage unification scenario. 



does not lie in an overlap region, the correct scenario would be selected over the 
other two. 

It is possible that other measurements (for example measuring Higgs couplings) 
will constrain some of the parameters of the model. For this reason, we will now 
assume that two of the four variables 9, </>, tan (3 or m 3 /2 are constrained in order to 
see the resulting improvement in string scenario discrimination. An obvious example 
to examine is the dilaton domination scenario (corresponding to 9 = 90°), which has 
attracted attention in the past in its own right. In figure 0a, we plot the ratios of 
masses x = m^ h jm^ R against y = m g /m u ~ L in the dilaton domination limit for all 
scanned values of the gravitino mass and tan (3 in table [1]. We can clearly distinguish 
between the three scenarios because there is no overlap in their regions in ratio 
parameter space. From the figure, we note that at least a 5 3% measurement of 
m g /m u ~ L and a 4% measurement in mg L /mg R would be required to always separate 
the models on the basis of these ratios alone. This is a very positive result: assuming 
dilaton domination, we have complete discrimination. Figure |7]b shows that there 
are minimal changes if one averages over the first two generations of squarks. The 
discriminating accuracy becomes smaller: 1% and 3% in rrig/mq avg and rrie L /me R 
respectively. 

Figure |8| illustrates departures from dilaton domination. m 3 / 2 and tan/3 have 
been fixed in the figure, whereas 9 and have been scanned over as in table |1|. The 
figure shows the variation of the ratios me R /m g and mq avg /m g as well as the variation 
of m H o/m~ b2 and m^o/m^ with string scenario. Figures [8]a and |8|b have tan/3 = 30 
and m 3 / 2 = 250 GeV fixed to different values than figures |c and |8|d (tan/? = 6 and 

5 We always quote the fractional uncertainty, e.g. on a ratio R, we quote the error required upon 
AR/R. 
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Figure 8: Highly discriminating ratios of sparticle masses for (a), (b): tan/3 = 30, = 
250 GeV and (c), (d): tan/? = 6, m 3/ / 2 = 200 GeV. Each plotted point corresponds to 
one (9, </>)-pair in the scan not excluded by the bounds. Green crosses are predicted by the 
GUT unification scale, early unification is denoted by red crosses, and blue stars are valid 
for the mirage unification scenario. 



m 3/2 — 200 GeV), in order to illustrate how the required discrimination accuracy 
can depend upon the supposedly constrained parameters. Discrimination is shown 
to generally be achieved for errors smaller than 14% in me R /m g and 4% in rriq avg /mg 
from figure ^]a. These numbers happen to be roughly identical for a different choice 
of tan (3 and m 3 / 2 , as shown in figure |8]c. Figures |8|b,d show how discrimination could 
be achieved with 5,7% errors respectively on m#o/m^ 2 and m^o/m^. We note that 
the overall values of these two ratios in any one scenario depends sensitively upon 
the values of tan/? and m 3 / 2 . tan/5 and m.3/2 could then be constrained by their 
measurement. 

We now fix the goldstino angles, to see if it is possible to distinguish the models 
when they are scanned over 777,3/2 and tan/5. The usefulness of combinations of 
ratios of sparticle masses depends in general upon the value of these fixed goldstino 
angles. As figures |^a-|9|b show, however, that the required accuracy on the ratios 
m qavg/ m g an d m e R / m g does vary with respect to a variation of the goldstino angles. 
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Figure 9: Discriminating ratios of sparticle masses in departures from dilaton domination: 
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(m 3 / 2 ,tan/3) pair in the scan not excluded by the limits. Green crosses are predicted by 
the GUT unification scale, early unification is denoted by red crosses, and blue stars are 
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From figures |9]a,b we see that for = 0°, {6 = 30°, 6 = 60°) we require a (3%, 5%) 
measurement of mq avg /rrig and a (11%, 10%) measurement of me R /mg to discriminate 
the three scenarios. In fact, these ratios discriminate for any fixed value of 9, <p but 
do not completely discriminate once the goldstino angles are scanned over. It is 
therefore difficult to say much that is definitive about the ratios in this case. It is 
possible that specific models will predict the goldstino angles, and in that case the two 
ratios mentioned above will be useful. However, the connection between measured 
quantities and 9, <fi is perhaps less obvious than for m 3 / 2 and tan/3. This means that 
it is difficult to see how one would directly infer their values from experimental data, 
and so we leave the discussion of constrained goldstino angles here. 

3.6 Interpretation 

Note that the required accuracy quoted here is the one necessary to separate the 
closest two points between different scenarios, i.e. assuming that the data lie within 
one of the predicted regions, what accuracy would definitely exclude the other two 
regions. This is the most pessimistic scenario, and we may well require less accuracy 
depending upon where the data finally lie. For example figure shows that the 
data m^ h /me R ~ 3.1 ± 0.6 and m g /mu L ~ 1.05 ± 0.03, would imply that the GUT- 
scale scenario would be selected over the other two whereas for m^/m^ ~ 1.8 and 
m g/ m a L ~ 1-0, uncertainties of (0.2,0.02) would be required respectively. Where 
the data finally lie will also decide the most important measurement. We take the 
dilaton dominated scenario to exemplify this point (see figure 0a). If the data lie 
over the GUT-scale region, the most important errors are those upon m^ L /me R as 
this will best separate the GUT-scale scenario from the other two. If, on the other 
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hand, the data lie over the early unification scenario, it is more important to have 
smaller uncertainties upon mg/mQ L for discrimination against the mirage unification 
scenario. 

We must ask what effect theoretical errors have on the interpretation of the 
results. Comparisons between the spectra derived from SOFTSUSY code and 
ISASUGRA [TS|], while showing qualitative agreement, display differences of order of 



the mass differences required to discriminate the scenarios |2g, ^5|. The errors are 
both dependent upon the flavour of particle (sleptons tend to have very good agree- 
ment, whereas coloured objects typically show a 3-5% difference) and on the region 
of parameter space. For extreme places in parameter space, for example near the ra- 
diative electroweak symmetry breaking boundary, predicted mass differences can be 
large between the two codes |^6, E5j . However, provided the point in parameter space 
being examined is not particularly special, any theoretical errors in codes should not 
change our conclusions about discrimination. This is because the errors due to (for 
example) not including finite threshold corrections in ISASUGRA will affect all points 
in the three scenarios in roughly the same way, moving them all in one direction for 
example. The statistical errors m t = 175 ± 5 GeV similarly provide an uncertainty. 
We note that this error should be drastically reduced by Run II of the Tevatron. 

It is clear, however, that better accuracy in the mass predictions will be necessary 
when interpreting real data. Of course, the discrimination depends at face value on 
the assumption of prior knowledge either of the goldstino angles (which will probably 
not be owned) or of m 3 / 2 and tan/3. The latter pair of parameters may be predicted 
in the hypothesis of each scenario, since gaugino masses can help us divine r/7.3/2 



and tan/? may well be obtained from Higgs measurements |24[. We note that each 
discriminating plot is valid for two particular values of fixed parameters. If the fixed 
parameters are not determined beforehand, they would have to be fit to the sparticle 
masses. This fitting procedure will be subject more to the theoretical errors, and 
will require them to be shrunk from their current levels. We would of course use all 
available (i.e. measured) highly discriminating ratios simultaneously in order to try 
to pin down the correct scenario. This would provide additional confirmation that 
the identified scenario is the correct one. 

It is likely that the next machine capable of producing TeV-scale superparticles 
in sufficient numbers to make a reasonable measurement is the LHC It remains to be 
seen how good the LHC is in determining ratios of masses, but it cannot constrain any 
single sparticle mass to a very accurate level (10% or so on the absolute values). The 
LHC can, however, constrain certain functions of masses to about 1% [^7], |28|, |7|, ^9f 



by finding the end-points of kinematic distributions. In ref. pSfl , the three most 
accurate measurements were found to be on the edges in the invariant masses detailed 
in table 0. They follow parts of the decay chain of q — > qx®, X2 ~ * ^> then I — > 
The existence of this chain depends upon the mass ordering rriq > m x o > mj. 



- 18 - 




400 

z max x 2 // maxx2 
(m llq ) /(m u ) 



1000 



Figure 10: Ratios of LHC edge variables. Each plotted point corresponds to one data 
point in four-dimensional parameter space (8, </>, tan/3, "13/2) that is not excluded by the 
imposed limits and where m x o > fn\ R - 



The edges shown in table || measure p9 



tot 



{miq(law) j 



(m 



llq 



(m 2 - m~ )(to~ 



I — toV 

xv 



min 



max 



to- 

(m 2 ~ — m 2 )(m 2 — to- ) (m^~m 2 () )(m 2 —m 2 ) 

K 1 xV y X 2 Ir' v 1 X2 JK Ir XV 



m 2 

(m~ — m 2 o)(m 2 o 

y " X2> K X2 



[valval 



m 2 

X°2 

m 2 m 2 )(m 2 — m 

Xl X2 JK X2 



xV 



1 1 

m~ TO n 
Ir X2 



2m l ~ m \\ 

(m 2 g - m 2 R ){m 2 R 



mi 



v 



(3.6) 



except for the special case in which to- 



when 



llq 



to„ 



m 



< m 2 -m 2 < m 4 and to 4 to 2 < m~m~ 

1 Xi X 2 Xi Xi 1 Ir 

fi) 2 instead. We use m^ ax to normalise the other edge 



variables, as its small error is negligible. 



We therefore plot the ratios (rnf^) 2 / (m^ ax ) 2 against (TOg(iow)) 2 /(^/; ax ) 2 m fig- 



ure 



10 in order to see if these precision LHC measurements can be used to distin- 



guish the string scenarios. All four free parameters are scanned over and all current 
experimental constraints (including g — 2 of the muon) are applied. Additionally, 
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we require m Y o > mr , to make sure that the above mentioned decay chain actually 



exists. This is the case in roughly one third of the investigated scan points. 

As figure [10] shows, there is no clear separation 
between all scenarios on the basis of these LHC edge 
variables. However, the (m^ (low) ) 2 /(m™) 2 -axis pro- 
vides some exclusion limits: the GUT scale case could 



edge variable 


accuracy 


m max 


0.1% 
1% 
1% 



Table 2: Accuracy of most 
useful LHC edge variables. 



be ruled out if (m; g(low) ) 2 /(m^ ax ) 2 does not lie between 
7 and 12. The horizontal axis appears to contain no 
discriminatory power. 



4. Conclusions 

We have seen that general information can be obtained about string theory scenarios. 
We considered three general string scenarios: two with the MSSM superfield content 
(GUT scale unification and mirage unification at the intermediate scale) and one with 
additional leptons and unification at the intermediate scale. We have mapped out 
the available parameter space for each scenario, using current empirical constraints. 

Our main results are better summarised in figures |] and || Figure |3| provides 
a map of viable parameter space slices for each of the three scenarios. We can 
see clearly that for different combinations of values of the goldstino angles different 
regions of 1713/2 — tan/3 parameter space are allowed. 

In figure ^| we find that ratios of particle masses may be able to separate each 
scenario when all parameters: 9, (f), 777.3/2, tan (3 are considered. This allows us to 
differentiate between the three scenarios provided the data do not appear in certain 
overlap regions. Clearly, fixing some of the parameters makes the separation between 
scenarios cleaner (see figs. 0-0) • We then determined what accuracy is required on 
the measurement of these masses in order to distinguish the scenarios. The errors 
are required to be less than the few percent level, depending upon the exact ratio. 

A previous study |30| found that the exotic heavy leptons present in the early 



unification scenario, can be discovered by the LHC if their masses are less than 980 
GeV. This could confirm the mirage scenario, making further discrimination redun- 
dant between the models considered here. The exotic leptons would not be discovered 
if they were heavier than 1 TeV, so all three scenarios would still be left to be dis- 
tinguished in that case. We note that a previous study confronted a fundamental 
SUGRA model with detailed linear collider and LHC measurements in a bottom-up 
approach Presumably, this kind of analysis could be repeated for stringy sce- 



narios and would complement the present study. It could be used to provide accurate 
measurements of high-scale soft-breaking parameters and to determine at what scale 
they unify. 

The small percent-level errors required on mass ratios are not expected to occur 
at the Tevatron experiments or even at the LHC M. But it is likely that a 
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Table 3: Lower bounds on relevant sparticle masses coming from the parameter scans, 
in GeV. The rows are marked by string scenario: early unification, mirage and GUT-scale 
unification respectively. 



future linear e + e~ collider facility [|4| would have sufficient accuracy, provided it 
had sufficient centre of mass energy to produce some of the particles involved in the 
ratios. For this reason, we present the lower bounds on relevant sparticle masses 
in table |3] for each scenario. Notice that this is an explicit way to differentiate 
the different scenarios. If, for instance, the neutralinos are discovered with a mass 
smaller than 60 GeV this would clearly favour the GUT scenario over the other two. 
Although the LHC can cover the lower bounds on each sparticle mass, only a multi 
TeV facility such as CLIC would provide enough centre of mass energy to produce 
sparticles if their masses are several TeV. We note, however, that lower values of 777.3/2 
(and therefore lower sparticle masses) are favoured by the fine-tuning parameter ||. 

There are several directions in which our work can be extended. First we may 
incorporate the charge and colour breaking (CCB) constraints on all the scenarios, 
as was done for the dilaton domination scenario in ref . || . While it would certainly 
be interesting to know if the universe were in some meta-stable vacuum, a global 
CCB vacuum might not necessarily rule the model out [p3]| . We could impose some 
arbitrary fine-tuning constraint on the models in order to restrict the maximum 
values of 7773/2; this would presumably increase the discriminatory power of our mass 
ratio tests. Another extension relates to the consideration of moduli domination, a 
limit that we did not consider since in that case the soft breaking terms appear only 
at the loop level and therefore anomaly mediation cannot be neglected. It is an open 
question to correctly combine anomaly mediation with gravity mediation for all soft 
breaking terms. 

This study is a crude first step towards the experimental discrimination of string 
models, leaving plenty of room for more investigation. For example, we only con- 
sidered two mass ratios at a time for discrimination, but this could be generalised 
to higher numbers, requiring an automated algorithm for finding the most discrim- 
inating ratios. It would be interesting to know the search reach in parameter space 



of the various future colliders as could be performed by HERWIG [34]. Further, we 



would like to know under which conditions experiments will be able to measure the 
discriminatory mass ratios. 

We hope that the techniques introduced in this article can be used to discrimi- 
nate other SUSY breaking scenarios, assuming eventual discovery of supersymmetric 
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particles. 
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